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Abstract 

We study the off-forward quark distributions (OFQD's) in the chiral quark-soliton 
model of the nucleon. This model is based on the large- iV c picture of the nucleon as 
a soliton of the effective chiral lagrangian and allows to calculate the leading twist 
quark- and antiquark distributions at a low normalization point. We demonstrate 
the consistency of the approach by checking various sum rules for the OFQD's. We 
present numerical estimates of the isosinglet distribution H(x,£, A 2 ). In contrast 
to other approaches we find a strong qualitative dependence on the longitudinal 
momentum transfer, £. In particular, H(x,£, A 2 ) as a function of x exhibits fast 
crossovers at \x\ = £/2. Such behaviour could lead to a considerable enhancement 
of the deeply-virtual Compton scattering cross section. 
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1 Introduction 



Although a large portion of quantitative information about strong interactions is con- 
tained in the well-known parton distribution functions of the nucleon, one should not 
forget that they provide us with still far from complete knowledge of the structure of the 
nucleon. Recently, a new type of parton distributions has attracted considerable interest, 
the so called off-forward parton distributions (OFPD's), which are generalizations simul- 
taneously of the usual parton distributions and of the elastic nucleon form factors. Taking 
the n-th moment of the OFPD's one obtains the form factors (i.e., non-forward matrix 
elements) of the spin-n, twist-two quark and gluon operators. On the other hand, in the 
forward limit the OFPD's reduce to the usual quark, antiquark and gluon distributions. 
In other words, the OFPD's interpolate between the traditional inclusive (parton distri- 
butions) and exclusive (form factors) characteristics of the nucleon and thus provide us 
with a considerable new amount of information on nucleon structure. 

The OFPD's are not accessible in standard inclusive measurements. They can, how- 
ever, be measured in deeply-virtual Compton scattering (DVCS) [1], 0, [| || |7j and in 
hard exclusive electroproduction of mesons 0, [|, A firm theoretical basis for the QCD 
analysis of DVCS and hard exclusive electroproduction of mesons in terms of off-forward 
parton distributions in the nucleon is provided by the recently proven factorization theo- 
rems for these processes ||, ||, 0] . A quantitative description of these classes of processes 
requires not only knowledge of the perturbative evolution of the OFPD's 0, |], [7], ||, but 
also non-perturbative information in the form of the OFPD's at some initial normalization 
point. Thus, a computation of the OFPD's at a low normalization point in a realistic 
model of the nucleon is of great importance. While the usual parton distributions have 
been measured in a variety of different experiments and can be confronted with model 
calculations a posteriori, in the case of the off-forward distributions the situation is rather 
opposite: Here, model calculations of the distributions at a low normalization point are 
required to determine the very feasibility of measuring the OFPD's, say, in DVCS at 
experimental energies. 

Recently it was shown that the chiral quark-soliton model of the nucleon [10], which is 
based on the instanton model of the QCD vacuum, provides a framework for a succesful 
calculation of the nucleon parton distributions, both unpolarized and polarized, at a 



low normalization point [IT], [L2], [13], [14]. This is a quantum field-theoretical description 
of the nucleon, with explicit quark degrees of freedom, which allows an unambiguous 
identification of the quark as well as antiquark distributions in the nucleon. We were 
able to demonstrate that all general properties of the quark and antiquark distributions 
(positivity, sum rules, Soffer inequalities, etc.) are correctly realized in this description. 
At the same time, the chiral quark-soliton model gives a good account of a variety of static 
properties of the nucleon and other baryons, and of baryon form factors (for a review, see 
ref. ||). Given the success in describing both nucleon form factors and the usual parton 
distribution functions, this approach seems well posed for a theoretical investigation of 
the off-forward parton distributions of the nucleon at a low normalization point. 

Recently, Ji et al. have studied the off-forward quark distributions in the bag model 
||, providing the first model estimates of these quantities. However, bag models encounter 
severe problems in applications to parton distribution functions. It is usually assumed 
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that the three quarks in the bag give rise only to quark distributions; however, they 
inevitably produce also an antiquark distribution with negative sign. This circumstance is 
not often emphasized. To overcome this problem one would need to add the contribution 
to the structure function arising from the forces which bind the quarks — in the case 
of the bag model, from the bag surface. Unfortunately, the latter is not described in 
terms of fields. In contrast to the bag model our approach ensures the positivity of 
both quark and antiquark distributions, because the forces which bind the quarks in 
the nucleon are described consistently in terms of fields. Below we shall see that in the 
case of off-forward quark distributions (OFQD's) the field-theoretic description leads to 
characteristic crossovers of the OFQD's at x = ±£/2, as opposed to the bag model result, 
which is smooth at these points. 

In this paper we study the off-forward quark and antiquark distributions in the chiral 
quark-soliton model of the nucleon. In section ^| we outline the QCD definition of the 
OFPD's and their basic properties. Section [3] gives a brief introduction to the chiral quark 
soliton model of the nucleon. In section ^ we derive the expressions for the OFQD's in the 
chiral quark soliton model and show that the sum rules are satisfied within the model. We 
discuss the contributions to the OFQD's from the bound-state level of quarks as well as 
from the Dirac continuum, whose presence is a consequence of the field-theoretic character 
of this description of the nucleon. We show that the latter is crucial near x = ±£/2 and 
lead to sharp crossovers of the OFQD's. Numerical estimates for the singlet unpolarized 
distribution H(x, £, A 2 ) are presented in section |5|. A summary and conclusions are given 
in section ||. 

2 QCD definition of off-forward quark distributions 

In QCD the off-forward parton distributions are defined through nondiagonal matrix 
elements of product of quark fields at light-cone separation. Here and in the following, 
we shall use the notations of refs. [1], |[]: 



dX 
— ( 

2vr 



-e 



■ iXx (P'\^j{-Xn/2)i^{Xn/2)\P) = H(x, £, A 2 ) U(P') iji U (P) 



+ ^E(x,^A 2 )U(P') i a^n,A u U(P). 



(2.1) 



Here n M is a light-cone vector, 



n 2 = 0, n-(P + P') = 2, (2.2) 

A the four-momentum transfer, 

A = P' - P, (2.3) 

Mn denotes the nucleon mass, and U(P) is a standard Dirac spinor. The off-forward 
quark distributions, H(x, £, A 2 ) and E(x, £, A 2 ), are regarded as functions of the variable 
x, the square of the four-momentum transfer, A 2 , and its longitudinal component 

£ = -(n-A). (2.4) 
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In the forward case, P = P', both A and £ are zero, and the second term on the r.h.s. 
of eq. ( |2. 1]) disappears. The function H becomes the usual parton distribution function, 

#(x,£ = 0,A 2 = 0) = q(x). (2.5) 

On the other hand, taking the first moment of eq.( [2.i|) one reduces the operator on the 
l.h.s. to the local vector current. The dependence of H and E on £ disappears, and the 
functions reduce to the usual electric and magnetic form factors of the nucleon, 

f^dxHfa^A*) = F X (A 2 ), (2.6) 

J dxE{x,£,A 2 ) = F 2 (A 2 ). (2.7) 

Taking higher moments of the distribution functions one obtains the form factors of the 
twist -2, spin-n operators. 



3 Chiral quark— soliton model of the nucleon 



Recently a new approach to the calculation of quark distribution functions has been 
developed |12| within the context of the chiral quark-soliton model of the nucleon [pTD| . In 
this paper we apply this approach to the calculation of OFQD's. It is essentially based 
on the 1/N C expansion. Although in reality the number of colours N c = 3, the academic 
limit of large N c is known to be a useful guideline. At large N c the nucleon is heavy and 
can be viewed as a classical soliton of the pion field [115], |i~6|. In this paper we work with 



the effective chiral action given by the functional integral over quarks in the background 



pion field pi [EL 0: 



exp (iS ee [7T{x)}) = / DipDip exp [i / d 4 x ip(i@ - MW 5 )ip 



U 



exp [iir a (x)r a ] 



U 15 



exp [iTX 



1 



75 



U ] . (3.1 



2 2 

Here if) is the quark field, M the effective quark mass, which is due to the spontaneous 
breakdown of chiral symmetry (generally speaking, it is momentum dependent), and U 
is the SU(2) chiral pion field. The effective chiral action given by (|3.8|) is known to 
contain automatically the Wess-Zumino term and the four-derivative Gasser-Leutwyler 
terms, with correct coefficients. Equation ( |3.8|) has been derived from the instanton model 
of the QCD vacuum [|r5| , pp| , which provides a natural mechanism of chiral symmetry 
breaking and enables one to express the dynamical mass M and the ultraviolet cutoff 
intrinsic in ( |3.8| ) through the Aqcd parameter. The ultraviolet regularization of the 
effective theory is provided by the specific momentum dependence of the mass, M(p 2 ), 
which drops to zero for momenta of order of the inverse instanton size in the instanton 
vacuum, 1/p ~ 600 MeV. For simplicity we shall neglect this momentum dependence in 
the general discussion; it will be taken into account again in the theoretical analysis and 
in the numerical estimates later. 
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An immediate application of the effective chiral theory ( |3.8| ) is the quark-soliton model 
of baryons of ref. fllO| |, which is in the spirit of the earlier works pl~l , |22fl . According to 



this model nucleons can be viewed as N c "valence" quarks bound by a self-consistent pion 
field (the "soliton") whose energy coincides with the aggregate energy of the quarks of 
the negative-energy Dirac continuum. Similarly to the Skyrme model large N c is needed 
as a parameter to justify the use of the mean-field approximation; however, the 1/N C - 
corrections can be — and, in some cases, have been — computed |J. 

Let us remind the reader how the nucleon is described in the effective low-energy 
theory (|3.8|). Integrating out the quarks in (|3.8|) one finds the effective chiral action, 



S cS [7i a (x)] = -N c Sp\ogD(U), D{U) = id -H(U), (3.9) 

where H(U) is the one-particle Dirac hamiltonian, 

H(U) = -ij°j k d k + M 7 °f/ 75 , (3.10) 

and Sp . . . denotes the functional trace. For a given time-independent pion field U = 
exp(i7r a (x)r a ) one can determine the spectrum of the Dirac hamiltonian, 

H$ n = E n <$> n . (3.11) 

It contains the upper and lower Dirac continua (distorted by the presence of the external 
pion field), and, in principle, also discrete bound-state level(s), if the pion field is strong 
enough. If the pion field has unity winding number, there is exactly one bound-state level 
which travels all the way from the upper to the lower Dirac continuum as one increases 
the spatial size of the pion field from zero to infinity ||10" |. We denote the energy of the 



discrete level as E\ ev , —M < E\ cv < M. One has to occupy this level to get a non-zero 
baryon number state. Since the pion field is colour blind, one can put N c quarks on that 
level in the antisymmetric state in colour. 

The limit of large N c allows us to use the mean-field approximation to find the nucleon 
mass. To get the nucleon mass one has to add N c E\ ev and the energy of the pion field. 
Since the effective chiral lagrangian is given by the determinant (|3.9|) the energy of the pion 
field coincides exactly with the aggregate energy of the lower Dirac continuum, the free 
continuum subtracted. The self-consistent pion field is thus found from the minimization 
of the functional [[HJ 

M N = min N c I E lev [U] + £ (E n [U] - E®) \ . (3.12) 

[ E n <0 J 

From symmetry considerations one looks for the minimum in a hedgehog ansatz: 

UJx) = exp [nr a (x)r a ] = exp \in a r a P(r)] , r = |x|, n = -, (3.13) 

r 

where P(r) is called the profile of the soliton. 

The minimum of the energy ( [3 . 1 2| ) is degenerate with respect to translations of the 



soliton in space and to rotations of the soliton field in ordinary and isospin space. For 
the hedgehog field ( |3.13| ) the two rotations are equivalent. The projection on a nucleon 
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state with given spin (S3) and isospin (T 3 ) components is obtained by integrating over all 



spin-isospin rotations, R |1JJ [TO 



(S = T,S 3 ,T 3 \...\S = T, S 3 , T 3 > = JdR 4>* S ^ T {R) ■ ■ ■ 47t 3 (R) • (3-14) 

Here 4>s^{R) is the rotational wave function of the nucleon given by the Wigner finite- 
rotation matrix IIIBI, ITU : 



4^ 3 (R) = V2STT(-l) T ^D s _^ S3 (R). (3.15) 

Analogously, the projection on a nucleon state with given momentum P is obtained by 
integrating over all shifts, X, of the soliton, 

(P'|...| P) = J d 3 Xe'( p '- p »' x ... (3.16) 

4 Off-forward quark distributions in the chiral quark— 
soliton model 

We now turn to the calculation of the off-forward quark distributions in the chiral quark- 
soliton model. This description of the nucleon is based on the l/iV c -expansion. At large 



N c the nucleon is heavy — its mass is 0(N C ). Thus for the large-iV c nucleon eq. ( 2.1 
simplifies as follows: 

j —e^(P',S' 3 m-\n/2)4ij(\n/2)\P,S 3 ) = 26 s , S3 H(x,£, A 2 ) (4.17) 

- -^6 3 ^(a fc ) 5 ^ 3 E(x,e,A 2 ), 
M N 

where S3, S' 3 denote the projections of the nucleon spin. 

Before computing the quark distribution functions we must determine the parametric 
order in 1/N C of the kinematical variables involved. Generally, when describing parton 
distributions in the large- N c limit, x ~ l/iV c , since the nucleon momentum is distributed 
among N c quarks. Furthermore, as in the calculation of nucleon form factors we consider 
momentum transfers A 2 ~ N®; hence, in particular, £ ~ l/iV c , so that £ is of the same 
parametric order 

The calculation of the off-forward parton distributions proceeds in much the same way 
as that of the usual parton distributions [T^]. One expands the quark fields in eq.( [4.17| ) in 



the basis of quark single-particle wave functions in the background pion field, eq. (|3.11| ). 
The nucleon matrix element is then obtained by summing over all occupied single-particle 
states (the bound-state level and the negative-energy continuum) and projecting on nu- 
cleon states with definite spin and momenta by integrating over collective coordinates of 
the soliton field with appropriate wave functions, cf. eqs.( [3.l4l , |3.16|) . When integrating 
over soliton rotations one must consider separately flavour singlet and nonsinglet matrix 
elements. In the flavour singlet case the integral over soliton rotations is trivial, and one 
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obtains 



J2 = S'X,S' 3 \ ^ / (-An/2)# / (An/2) |P,T = S,T 3 ,S 3 > 

f 

Wss'tW^ / d 3 X exp[iA ■ X] 



iV,M 2 



7T 

x X! exp{z^ [(x + ^/2)Af w -E n ]}^(X)7 #n(X-z e3), (4.18) 

occup. 



where e 3 is the unit vector in the third direction. This equation was derived in |T2[ for the 
case P' = P; the generalization to P' ^ P is straightforward. [Here and in the analogous 
expressions below it is understood that one subtracts the sums over levels of the vacuum 
hamiltonian (U = 1).] For the flavour nonsinglet part, on the other hand, one has 

d\ elXx r = g,^ # | ^(_ An /2)r a #(An/2) \P,T = S,T 3 ,S 3 ) 
2n 



N C M 2 N 



t 'dz° J d 3 X exp[*A ■ X] J dR(f)^f*{R)(f^= s s 3 {R) -Tr^VPr 6 ) 
x Yl exp{i^[(ar + e/2)M JV -S n ]}$t(x)r 6 7V$ n (X-2 e 3 ). (4.19) 

occup. 



Here <Pt s s^(^ l ) are the rotational wave functions of the soliton, eq.( |3.15 ). The integral 
over the soliton orientation matrix R can be computed, 



dR^l s ,(R)<f>l S3 (R) -Tr(rfr a RT b 



( T )t't 3 {o- )s's 3 



(4.20) 



Comparing eq. (|4.18| ) and eqs. (|4.19| , |4.20|) with eq.( |4.17| ) we immediately see that in the 
leading order of the 1/N C expansion only the flavour singlet part of if (x, £, A 2 ) and the 
flavour-nonsinglet part of E(x,£, A 2 ) are non-zero. They are given, respectively, by 



N C M N 



2vr 



J dz° J d 3 X exp[iA ■ X] 



E^/(^e,A 2 ) 

/ 

x £ exp{^°[(x + e/2)M Ar - J E; ri ,]}<l>t t (X)(l+7V) < f n .(X-^e3), (4.21) 

occup. 

e 3i*A^ 3 W,A 2 ) = - l J^Ml [ dz o [ rf 3 X exp r, A . X ] 

3n J J 

x J2 exp{tz [(x + ^/2)M N -E n }}^ n CX)r k (l + 1 1 3 )^ n CX-z e 3 ). (4.22) 

occup. 

The isovector part of H(x, £,A 2 ) and the isosinglet part of E(x, £, A 2 ) appear only in 
the next-to-leading order of the l/iV c -expansion, i.e., after taking into account the finite 
angular velocity of the soliton rotation. 

Before going ahead with the evaluation of the expressions eqs.( [OT , |4.22| ) we would 
like to demonstrate that the two limiting cases of the off-forward distributions — usual 
parton distributions and elastic form factors — are correctly reproduced within the chiral 
quark-soliton model. Taking in eq.( |4.21| ) the forward limit, A — ► 0, one recovers the 



formula for the usual singlet (anti-) quark distributions in our model which was obtained 
in ref.[|l2|- Thus the forward limit, eq. Q2.5Q , is reproduced. On the other hand, integrating 
eqs.( [4.21| , [4.22D over — 1 < x < 1 one obtains (up to corrections parametrically small 
in 1/N C ) the expressions for the electromagnetic formfactors of the nucleon derived in 
ref. 



J^dxY.Hfix^^ 2 ) = N c Jd 3 Xexp[iA-X] £ $+(X)$ n (X) 



occup. 



Ff =0) (A 2 ), (4.23) 
iN c M N 



j_teE®{x,t-,£?) = --^^e^Jd'Xex^A.X] £ ^(XhV^X) 

occup. 

= F 2 (T=1) (A 2 ). (4.24) 



The electromagnetic formfactors computed in the chiral quark soliton model on the basis 
of these formulas compare very well with the experimentally measured ones up to momenta 
of order A 2 ~ 1 GeV 2 |,|23J. 

Eqs. Q4.21|J4.22| ) express the OFPD's as a sum over quark single-particle levels in the 



soliton field. This sum runs over all occupied levels, including both the discrete bound- 
state level and the negative Dirac continuum. We remind the reader that in the case of 
usual parton distributions it was demonstrated that in order to ensure the positivity of 
the antiquark distributions it is essential to take into account the contributions of all oc- 
cupied levels of the Dirac Hamiltonian |L2). The so-called "valence level aproximation" for 



structure functions in the chiral quark-soliton model advocated in leads to unaccept- 
able negative antiquark distributions. We shall see below that also in the off-forward case 
the contribution of the Dirac continuum drastically changes the shape of the distribution 
function, leading to characteristic crossovers of H(x,£, A) at |x| = £/2. 

We shall now compute the contributions of the discrete bound-state level and the 
negative Dirac continuum to eqs. ( [4.21 , 02 ). We focus here on the isosinglet distribution 



H(x,£,A); the discussion and the expressions derived below can be easily generalized to 
the case of other OFQD's. 

The contribution of the discrete bound-state level to eq. Q4.21Q can be computed using 



the expressions given in the Appendix. The result is shown in Fig. 1 for the forward case 
and Fig. 2 for a non-zero momentum transfer. Being taken by itself this contribution 
resembles the shape of H(x, £, A 2 ) obtained in the bag model p|. 

To calculate the contribution of the Dirac continuum to eq. ( [4.21 ) we resort to an 



approximation which proved to be very successful in the computation of usual parton 



distributions, the so-called interpolation formula [12]. One first expresses the continuum 
contribution as a functional trace involving the quark propagator in the background pion 
field. The quark propagator can then be expanded in powers of the formal parameter 
dU/(—d 2 + M 2 ), which becomes small in three limiting cases: i) low momenta, \dU\ <C M, 
ii) high momenta, \dU\ 3> M, Hi) any momenta but small pion fields, | log f/| <C 1. One 
may therefore expect that this approximation has good accuracy also in the general case. 
As was shown in refs. [ITSL |TJ] for usual parton distributions this approximation preserves 
the positivity of the antiquark distributions and all sum rules; moreover, it gives results 
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very close to those obtained by exact numerical diagonalization of the Dirac hamiltonian 
and summation over the negative-energy levels. 

To derive the interpolation formula for the isosinglet distribution H (x, £, A 2 ) (the ex- 
pressions can easily be generalized to the case of isovector E^ 3 '(x, £, A 2 )) we proceed in 
analogy to the case of usual parton distributions [I2| and rewrite the continuum contri- 
bution to eq. Q4.21Q as 



£#/(*,£, A 2 ) 
/ 



N C M N 
2nT 



Im / dz e 



Az°xM N 



cont 



xSp (v^) exp[*A • X + z\v^)] 



1 



[i^d^ - MUt* + iO] 



(tf -> 1) . (4.25) 



Here we have introduced a dimensionless light-like vector, = (1,0,0, 1), and T is the 
time interval which is canceled by a corresponding factor arising from the functional 
trace. We have written here explicitly the vacuum subtraction term. Expanding now in 
dU / {—d 2 + M 2 ) and evaluating the the functional trace in momentum space one obtains 
in leading order: 



A 2 ) 



cont 



-2M N N C Im 
M(p 2 ) 



d 3 k r d 4 p 



5 



x 



(2tt) 3 7 (2tt) 4 
M((p-k) 2 



(p 2 - M 2 + iO) (p - k) 2 - Ml + iO 

e^-e, a- -a), 



(x - -)M N - v ■ p 

(k-v)Tr a .{u(k-A)[U(k)] + } 



(4.26) 



where the Fourier transform of the soliton field is defined as 



£/(k) = / ^xe- !k x [E7(x) - 1] 



(4.27) 



In eq. ( [4.26] ) we have re-instated the momentum dependence of the constituent quark 
mass, M(p 2 ), which cuts the loop momentum p and thus regularizes the UV divergence; 
by M we denote the value of the mass at zero momentum, M(0). We have neglected in 
eq. ( |4.26| ) the momentum dependence of the quark masses appearing in the denominators; 
the masses standing in the denominators do not play the role of an UV regulator and their 
momentum dependence is not essential. 

In ref. ||12|| the continuum contribution to the quark distribution functions was com- 
puted by regularizing the loop integrals with a relativistic Pauli-Villars cutoff. One may 
argue that this regularization mimics the effect of the momentum dependence of the con- 
stituent quark mass. [We shall soon see under which conditions this assumption is valid.] 
Let us evaluate the continuum contribution to the OFQD, eq. Q4.26 ), also with a Pauli- 
Villars cutoff, neglecting the momentum dependence of the masses in the numerator. One 
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obtains 



A 2 

. / 

d 3 k 



N c MnM% 



2tt 



+ K 




sign x 



x-^\M N -¥ 



ReTr /z {C/(k-A) U(k) 



2 (x - |) Mjv [(a; - §) - P 



(4.28) 



where the last integral over transverse momenta is supposed to be regularized by the 
Pauli-Villars method, i.e., by subtracting the equivalent integral with M replaced by the 
regulator mass, Mpy; see for details. [Here 6 denotes the step function.] Eq. ( |4.28| ) has 
the remarkable property of being discontinuous at |x| = £/2. Such behaviour would lead 
immediately to a violation of the factorization theorems for, say, deeply virtual Compton 
scattering (DVCS) processes^ 0, since the expression for the DVCS amplitude contains 
a factor 



Re f 1 dx \ y + ) ) £#(x,£, A 2 ), (4.29) 

J-i \x-£/2 + i0 x + Z/2 + iOj 



which would be logarithmically divergent if H(x, £, A 2 ) were discontinuous at \x\ = £/2. 
However, this conclusion is premature: the discontinuities of eq. Q4.28| ) are artifacts of 



neglecting the momentum dependence of the constituent quark mass, as we shall now 
show. 

Let us analyze the original loop integral over p in eq. ( [4.261 ). By inspection of the 



denominators it is easy to see that for x close to ±£/2 the dominant contribution to 
the integral comes from the region where the virtuality in one of the quark propagators 
becomes large, 

However, in this region, due to the presence of the momentum-dependent mass M((p—k) 2 ) 
in the numerator, the whole expression drops to zero. More precisely, one can see that for 

M-f/2 ~ (4.31) 

M N 

where p ~ (600 MeV) -1 is the characteristic scale of momentum dependence of the quark 
mass (the average instanton size), the integral over p± is cut already at p\ ~ Mq, whereas 

2 We are grateful to A.V. Radyushkin for discussion of this point 
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with a momentum-independent mass it would always be cut at p\ ~ 1/p 2 , see eq. ( |4.28| ). 
This analysis shows that the momentum dependence of the constituent mass can not be 
neglected for x in region eq. ( |4.31| ) . Consequently, keeping the momentum dependence 
of the constituent quark mass, the discontinuities present in eq. (|4.28| ) are smeared over 
an interval of order (M^p) 2 /N c , so that the integral ( |4.29| ) is finite but it is still para- 
metrically (and numerically) enhanced by a factor of log[(Mop) 2 /iV c ]; see the next section 
for numerical estimates. For values of x far from ±£/2 the momentum dependence of 
the mass can be safely neglected and the simplified expression eq. ( |4 . 2 8| ) gives reliable a 
approximation to eq.( [4.26| ). 

Let us also note that the points x = ±£/2 divide the interval of the variable x (— 1 < 
x < 1) in three regions: x < — £/2, where the function H(x, £, A 2 ) describes the antiquark 
distribution; x > £/2, where it corresponds to the quark distribution, and — £ < x < £/2, 
where H(x, £, A 2 ) resembles a meson wave function. It is therefore natural that the 
functions if(x, £, A 2 ) has crossovers at |x| = £/2. This feature is remarkably reproduced 
in our model owing to the Dirac continuum contribution. We note that in the forward 
limit, A — > 0, this crossover corresponds to the fact that both the quark and antiquark 
distributions are positive; hence the universal function q(x) (which describes the quark 
distribution at x > and minus the antiquark distribution at x < 0) evidently must 
exhibit a crossover at x — 0. 



5 Numerical results and discussion 

We have calculated numerically the isosinglet distribution H (x, £, A 2 ); the analogous cal- 
culations for other OFQD's can be done along the same lines. Our main purpose here 
is to discuss qualitative behaviour of the OFQD's; to this end it is sufficient to consider 
H(x,!;,A 2 ). A comprehensive study of the other distributions will be given elsewhere. 

For the numerical calculations we shall use the variational estimate of the soliton 
profile, eq.(eq. (EU3p), of ref.0 (M = 350 MeV), 

P(r) = -2 arctan ( ^| J , r w 1.0/M , M N w 1170 MeV, (5.32) 



which has been used in the calculation of usual parton distributions in refs. fll2|, p^j . Fur- 



thermore, we approximate the momentum-dependent mass predicted by the instanton 



model of the QCD vacuum [19] by the simple form 



where the parameter A is related to the averaged instanton size, p, by A = 6 1/3 p- 1 . This 
form reproduces the asymptotic behaviour of M(p 2 ) at large euclidean p 2 obtained in the 
instanton vacuum, 

M (~P ~ (P -> oo). 

p°p° 
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We have explored also other forms of the momentum dependence of the mass and found 
that numerically the results are very close to each other. 

We estimate the Dirac continuum contribution to if(x, £, A 2 ) using the interpolation 
formula, eq. ( |4.26|) , which gives a reliable approximation preserving all qualitative features 
of the continuum contribution. The contribution of the discrete level is calculated using 
eq. ([Op . 

First we compute H(x, £, A 2 ) in the forward limit, A —>■ 0, where it coincides with 
the usual quark and antiquark distributions. The result is shown in Fig. 1, where we plot 
separately the contributions of the discrete level and that of Dirac continuum (computed 
from the interpolation formula), as well as their sum. It should be emphasized that 
the distribution of antiquarks arising from the discrete level (see eq. ( |A.5|) ) is definitely 
negative and sizeable. []. Positivity of the parton distributions is restored only when 
one includes the contribution of the Dirac continuum. The full result for the function 
H(x, 0, 0) (consisting of the level and continuum contributions) exhibits strong crossover 
at x = 0, corresponding to the fact that both the quark and antiquark distributions are 
positive. The crossover occurs in an x-interval of order \x\ <~ 0.05; we shall see that 
this remains so also for the off-forward case. To illustrate this we show in Fig. 2 the 
distribution if(x,£, A 2 ) as a function of x, for £ = 0.3 and Af, = —A 2 — £ 2 \ l\ = 0. 
Again we have plotted separately the contributions of the discrete level and the Dirac 
continuum (according to the interpolation formula), as well as their sum. We see that 
the discrete level contribution is a smooth function which does not "know" about the 
points x = ±£/2. This kind of behaviour has been assumed in all model approaches to 
OFQD's, e.g. in the bag model calculations ||. However the contribution of the Dirac 
continuum changes drastically the picture: the function H(x, £, A 2 ) shows now a fast 
crossover around the points x = ±£/2. The interval of x over which this crossover occurs 
is of the order of ~ 0.05, as in the forward case (see Fig. 1). 

The fast crossover of H(x,£, A 2 ) at x — ±£/2. may have physical implications. For 
example, it may lead to the considerable enhancement of the DVCS amplitude, because 
the crossover occurs exactly at the points where the integral eq. (|4.29|) entering the DVCS 
amplitude has singularities. Numerically, at £ = 0.3 and A^ = the contribution of the 
Dirac continuum to the integral ( |4.29| ) is 5.1 (almost the whole integral is collected in 
small vicinity of point x = ±£/2), while the corresponding contribution of the discrete 
level is 3.2. We thus see that the crossover in H(x,£,A) contributes more than 60% to 
the integral ( 4.29| ) and may thus considerably increase the DVCS cross section. 

In order to illustrate the dependence of H (x, £, A 2 ) on £ and A 2 we plot this function 
for a fixed momentum transfer of A 2 = —0.5 GeV 2 for various values of £ (see Fig. 3), 
and for fixed £ = 0.3 and various values of momentum transfer (see Fig. 4). 

In the large N c limit the nucleon is heavy, so the OFPD's do not automatically go to 
zero at x = 1; see the discussion in [JT^] . However, at x ^> 1/N C the distributions behave 
as ~ exp(— const ■ N c x). Numerically, even for N c = 3 all distributions computed are very 
small at x « 1. 



3 In the extreme case of a very strongly bound discrete level when it approaches the lower continuum, 
this level would not produce quarks at all - only antiquarks, but with a negative sign! 
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6 Conclusions 



In this paper we have investigated the off-forward quark distribution (OFQD's) functions 
in the nucleon in the large- N c limit. At large N c the nucleon can be viewed as a heavy 
semiclassical body whose N c "valence" quarks are bound by a self-consistent pion field. 
The energy of the pion field is given by the effective chiral lagrangian and coincides with 
the aggregate energy of the Dirac sea of quarks (the free continuum subtracted). To 
compute the quark and antiquark distributions (forward and off-forward) one must sum 
the contributions from the discrete level and from the (distorted) negative-energy Dirac 
continuum. 

We have found that the flavour singlet off-forward distribution H(x, £, A 2 ) exhibits a 
qualitatively new behaviour due to the contribution of the Dirac continuum: The function 
shows fast crossovers around the points x = ±£/2. Our numerical estimates indicate that 
this crossover could lead to a considerable enhancement of the deeply virtual Compton 
scattering cross section. 

In the crossover regions the behavior of the OFQD's is essentially determined by the 
momentum dependence of the constituent quark mass generated in the dynamical breaking 
of chiral symmetry. In particular, the region in x over which the crossover takes place is 
proportional to (M p) 2 / N c , where p is the characteristic momentum scale at which the 
dynamical quark mass drops to zero — the average instanton size. Thus, the enhancement 
of the deeply virtual Compton scattering cross section is governed by the small parameter 
intrinsic to the instanton model of the QCD vacuum, the packing fraction of the instanton 
medium, (M p) 2 oc (p/R) 4 . 

The OFQD calculated here in the chiral quark-soliton model refers to a low normal- 
ization point of order the UV cutoff intrinsic in this model, ~ 600 MeV. To obtain the 
OFQD's at higher normalization points one has to evolve the "primordial" distributions 
using the evolution equations derived in refs.p, f|, [?]. With regard to the description 
of experiments it would be extremely interesting to study to what extent the crossovers 
at x = ±£/2 persist when evolved to higher normalization points. 
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A Bound-state level contribution to H(x, £, A 1 

We present here the contributions of the discrete bound-state level to the singlet H(x, £, A 2 ). 
The bound-state level occurs in the grand spin K — and parity II = + sector of the 
Dirac hamiltonian ( ft.lOj) . In that sector the eigenvalue equation takes the form: 



/ 8 2 

1 McosP(r) MsinP(r) 

or r 



d 

, M sin P(r) 

\ or 



h [r) 



-M cos P(r) 



J 



E\ 



lev 



h (r) 



(a.i; 



We assume that the radial wave functions are normalized by the condition 



drr 2 [h 2 (r) + jf(r)\ 



1. 



(A.2) 



We introduce the Fourier transforms of the radial wave functions, 



CXJ CXJ 

h(k) = J drr 2 h {r)R k0 {r), j{k) = j ' dr r 2 ji{r)R kl {r), (A.3) 



where 



RkAr) 



z / 2 r l ( 1 d\ sin kr 



(A.i) 



The bound-state level contribution to the singlet H(x, A 2 ) distribution function can 
be simply obtained from the general eq. ( |4.21| ). We get: 



where 



J2H f (x,^A 2 ) lev = 27rN c M N J 



d 2 ki 



f 

j(k)j(k 



{ 2,fkk\ mm 
„ k ± • (k ± + A ± ) + (xM N - E lcv ) 2 - \i 2 M 2 N 



- h(k)j(k' 



xM 



N 



kk 

E lev + \£M N 



k 



- h(k')j(k) 



xM N - E lev - UM. 



N 



k> 



k' 



'k ± 2 + ((x + -i)M N - E lcv y 

I (k ± + A T ) 2 + [(x - ^)M N - E, cv ] 2 



(A.5) 



(A.6) 



(A.7) 



Note that the r.h.s. of eq. (|A.5|) is positive in the forward limit (A — > 0) for all values of x, 
in particular at x < where eq. ( |A.5| ) determines in fact the antiquark distribution. Since 
q(x) = —H(—x,0,0), it means that eq. ( |A.5| ) gives a negative distribution of antiquarks 
at x > 0. At the same time it is easy to check by integrating over x eq. ( |A.5|) that the 
baryon number sum rule is fully saturated by the discrete-level contribution only. 
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Figure 1: The isosinglet distribution H(x, £, A 2 ) in the forward limit, A = 0. Dashed 
line: contribution from the discrete level. Dashed- dotted line: contribution from the Dirac 
continuum according to the interpolation formula, eq. (|4.26 ). Solid line: total distribution 
(sum of the dashed and dashed-dotted curves). 
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Figure 2: The isosinglet distribution H(x,^, A 2 ) for Ay = (A T = -A 2 - £ 2 M%) 
and C, = 0.3. Dashed line: contribution from the discrete level. Dashed- dotted line: 
contribution from the Dirac continuum according to the interpolation formula, eq. ( |4.26|) . 
Solid line: the total distribution (sum of the dashed and dashed-dotted curves). The 
vertical lines mark the crossover points x = ±£/2. 
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Figure 3: The isosinglet distribution H(x,C,,A 2 ) (total result) for fixed A 2 = -0.5 GeV 2 
and various values of £. 
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